求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

超音速燃烧二元流场的数值模拟

由二元超音速燃烧流场垂直喷射Ｈ２的数值模拟计算，描述了缝隙喷嘴附近氧气与空气的混合和燃烧过程，使用两步化学反应，ＭａｃＣｏｒｍａｒｋ显式时间分裂法和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ修正的代数湍流模型，计算超音速燃烧流场，得出了流线，静压，静温，Ｈ２和Ｈ２Ｏ浓度等值线图，反映了缝隙喷嘴附近回流区的大小，自动着火点的位置及火焰向主流的传播现象。     

降落伞充气过程的数值模拟

 基于降落伞的重要应用与设计的实际需要，降落伞的数值模拟开始得到越来越多的重视，而充气过程是其中最为复杂的一个阶段。本文建立了平面圆形伞主充气过程中的CFD（Computational Fluid Dynamics）与结构动力学的MSD（Mass Spring Damper）之间的耦合模型。流场求解采用稳定性较高的标准k-ε模型，在多块贴体坐标下，获得某时间节点处的流场，并将该流场中的压力数据引入MSD模型，以获得下一时间节点的伞衣形状，最终获得主充气过程中伞衣形状和流场之间的动态关系。数值计算结果和实验结果及经验值比较，均有较好的一致性。充气过程的数值求解有助于提高对降落伞充气过程机理的理解。      

液压能源管路系统振动主动控制的理论研究

对飞机液压能源管路系统的振动问题进行了探讨.提出一种基于振动主动控制(VAC)技术对飞机液压能源管路系统进行消振的方法,在手段上采用压电/压电致伸技术,以压电陶瓷(PZT)作为作动器,具有驱动力大、响应频率高和体积小等优点;控制方法上采用参数寻优的控制策略,它能始终能根据外界干扰因素的变化来调整控制参数.经过理论分析和仿真验证,证明其具有良好的自适应性和鲁棒性,在变转速和变负载情况下始终能保持最佳的消振效果.     

塞式喷管三维流场的数值模拟

从曲线坐标下的三维平均雷诺的N-S方程出发,并用k-ε两方程湍流模型封闭方程组,采用二阶精度的精度无波动、无自由参数的耗散差分格式(NND格式),发展了模拟塞式喷管三维流场的数值程序,并对线性塞式喷管和瓦状塞式喷管在不同落压比下的三维流场进行了数值模拟研究.结果较好地预示了塞式喷管的三维流场特征,并且显示瓦状塞式喷管除了圆柱形结构比平板更能承载外,其塞锥压强分布可以进一步改善其受力情况.     

半浮区液桥自由面振荡特征

使用一套双路同步显微放大录像及图像处理系统,研究了夹角为90°的两个竖直剖面沿轴向不同高度液桥自由表面的位移振荡,得到轴向不同高度上的环向相位差,观测到在振荡周期内液桥表面轮廓线在不同时刻均不相交的特点。实验结果表明,液桥自由面振荡除具有径向扰动及轴向振荡变化外,还存在环向波动变化。根据实验结果,讨论了环向波动基本变化规律以及本实验在测量精度和测量误差方面存在的问题。     

几种流体壁面剪切力测压管的性能研究

对四种壁画型流体剪切力测量差压计Ｓｔａｎｔｏｎ管，Ｋ型管，穴型管，及微孔表面差压管在湍流边界层条件下进行试验并对其结构特点进行比较后表明，这四种差压计具有相同的特征，所测量的区差均为流量动量变化所致，其流体剪切力可以表示为τ＝Ｃ△Ｐ￣ｎ，为避免对标定数据过于繁复的多段拟会，并使差压计能在较大测量范围内保持拟合公式的精度，各类差压管在与固体壁画铅垂方向上的尺度值应尽量低，微孔表面差压管结构复杂，灵敏度低，但它具有线性特性。     

粘性流体大幅晃动的ALE有限元模拟

采用有限元法数值求解了具有自由液面大幅移动边界的Navier－Stokes方程。对流体区域采用了任意的拉洛朗一欧拉（ALE）运动学描述，网格结点可以任意移动而不依赖于流体的运动，结合了拉格朗目描述易于处理移动边界与区域变形的优点和欧拉描述可克服单元缠结的优点，提出了简单合理的网格更新方法，能精确地跟踪运动的自由液面。为了更精确地处理强对流项的作用，采用了迎风流线Petrov－Galerkin（SUPG）加权余量法建立有限元方程。算例表明本方法对贮箱内流体的非稳态大幅晃动过程的数值模拟非常成功。     

热完全气体的热力学特性及其N-S方程的求解

首先用五次多项式拟合给出了温度在50～3000K范围内的热完全空气的焓值与温度之间的函数关系式，导出了其它热力参数e，cp，cv和γ的表达式。接着提出了热完全空气总温、总压的计算方法，并将其计算结果与量热完全空气的结果进行了比较。最后，将该热完全空气模型用于N-S方程求解，对NAPA软件进行了改进，并用该软件时高马赫数钝体绕流流场和乘波体流场进行了计算，分析了气体真实效应时流场结构及参数的影响。结果表明，本文提出的总温、总压计算方法及时NAPA软件的改进是成功的，可较准确地模拟高超声速流动的主要特征。     

建筑物风环境三维定常数值仿真

针对建筑物的风环境问题，采用计算流体动力学CFD软件对两幢建筑物进行了数值模拟，得出了建筑物的风场分布，以及风载的分布规律，结果表明，气流在来流方向第一幢建筑物的顶部、两侧的前沿产生分离，建筑物之间产生强烈的漩涡，顶部的局部负压最大。